Negative refraction and superlensing in a 2D photonic crystal structure 
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We experimentally and theoretically studied a new left-handed (LH) structure based on a photonic 
crystal (PC) with a negative refractive index. The structure consists of triangular array of rectan- 
gular dielectric bars with dielectric constant 9.61. Experimental and theoretical results demonstrate 
the negative refraction and the superlensing phenomena in the microwave regime. The results show 
high transmission for our structure for a wide range of incident angles. Furthermore, surface termi- 
nation within a specific cut of the structure excite surface waves at the interface between air and PC 
and allow the reconstruction of evanescent waves for a better focus and better transmission. The 
normalized average field intensity calculated in both the source and image planes shows almost the 
same full width at half maximum for the source and the focused beam. 



Left-Handed-Materials (LHM) are materials with si- 
multaneously negative dielectric permittivity e (r) and 
negative magnetic permeability fi(r). The phase veloc- 
ity of the light wave propagating inside this medium 
is pointed in the opposite direction of the energy flow. 
Thus, the Poynting vector and wave vector are antipar- 
allel, consequently, the light is refracted negatively. The 
existence and impact of such materials was first pointed 
out by Veselagc 1 . Years later, several theoretical and ex- 
perimental groups investigated the LHM 2 * 3 *^. For LHM 
based on PCs, Notomi studied light propagation in a 
strongly modulated two-dimensional (2D) PC 6 . Luo et al 
have studied subwavelength imagining in PC 7 . Cubukcu 
et al demonstrated experimentally single-beam negative 
refraction and superlensing in the valence band of 2D PC 
operating in the micro wave regime 8 . Foteinopoulou et al 
emphasized the time evolution of an EM wave as it hits 
the interface between a right-handed (RH) and a LH ma- 
terial interface 9 . Pendry was the first who suggested that 
a slab of lossless LHM with both e (r) and /i (r) equal to 
-1 should behave as a perfect lens 10 , i.e, the small de- 
tails as well as the larger ones could be reproduced by 
such a lens. The reconstruction of the evanescent wave 
components is the key to such perfect focusing. It is 
shown in another study 11 how the evanescent waves get 
amplified upon reaching the interface between a RH and 
LH medium and consequently how they participate in 
improving the quality of the image. Homogeneous meta- 
materials and PCs were used to demonstrate the effect. 
The main challenge is to find a structure with n = — 1 
for which the matching condition are verified with no 
reflections. The candidate structure should have three 
characteristics. First, it should exhibit almost isotropic 
Equal Frequency Surfaces (EFS) in the band with neg- 
ative curvature 5 . Secondly, it should guarantee the ab- 
sence of any higher order Bragg reflections for any angle 
of incidence. Finally it should guarantee single beam 
propagation. 

In this paper we propose such a structure based on 
PC. The proposed structure meets all the required con- 
ditions and has a refractive index very close to n = — 1 



as well as a high transmission at the desired frequency. 
Additionally, with an appropriate cut and a specific ter- 
mination of the surface, the proposed structure focuses 
the image in a better way by exciting the surface waves 
and enabling them to contribute toward the quality of 
the image. Furthermore, for different incident angles, 
the cut structure demonstrates a high transmission as 
well as a maximum coupling with the surface modes at 
the interfaces. The proposed structure consists of a tri- 
angular array of dielectric bars with a dielectric constant 
e = 9.61 in air. The dimensions of each bar in the x and 
y directions are respectively r x = 0.40a and r y = 0.80a 
where a = 1.5875 cm is the center to center separation 
between bars. The length of each bar is I = 45.72 cm. 
Fig. 1 shows the structure with more details. Only the 
transverse magnetic (TM) modes are considered here ( 
the E field is parallel to the rods). 

The plane wave method is used to compute the pho- 
tonic band structure as well as the EFS. Some of the 
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Figure 1: (a) A triangular array of dielectric bars in air with 
r x = 0.40a and r y = 0.80a, where a = 1.5875 cm is the lattice 
constant, (b) A picture of the structure with 33 layers in the 
lateral direction and 9 layers in the propagation direction. 
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theoretical results are obtained using the Finite Differ- 
ence Time Domain (FDTD) method. This method shows 
a time and space evolution of the emitted EM waves. 
More details of this algorithm can be found in reference^. 
In all the FDTD simulations we report, LI AO bound- 
ary conditions 13 are used. Depending on the desired 
conditions for the case under study, the source emits 
a monochromatic TM- (E-) polarization of desired di- 
mensionless frequency. A Gaussian source (Gaussian in 
space and "almost" monochromatic) is placed outside the 
structure to check the negative refraction. However, a 
point source also excites the surface waves at the inter- 
face between RH and LH material. The photonic band 
structure as well as the EFS in k space are shown in 
Fig. 2. The quantity (fa/c) is the dimensionless fre- 
quency where c is the velocity of light and a is the lattice 
constant. The second band, extending from 0.27 to 0.40, 
has a negative curvature and is of particular interest. To 
study the system and to compute the refractive index, 
the EFS are plotted for different frequencies. Note that 
the EFS consist of the allowed propagation modes for a 
specific frequency and in our case, they have closed con- 
cave like form. Within the second band the shape of the 
EFS tends to became circular for higher frequencies and 
reproduces the symmetry of the system for lower frequen- 
cies. It is important to notice that in contrary with the 
triangular lattice with circular cross section for which 
the hexagonal symmetry is conserved on the lattice as 
well as on the scatterers, this structure with rectangular 
rods lose its hexagonal symmetry for the scatters when 
the wave is traveling along a specific direction. Thus 
the breaking of the symmetry gives some anisotropy in 
the structure. Thus, the EFS are not totally isotropic 
because of the breaking of the symmetry. Therefore, the 
computed refractive index of our structure using the EFS 
plot is somehow affected by the anisotropy and is not ex- 
actly equal to —1 instead, it ranges from —0.97 to —1.22. 

Transmission measurements are performed to verify 
and test the negative refraction and superlensing. The 
experimental setup consists of an HP 85 10C network ana- 
lyzer, a waveguide horn antenna as the transmitter, and 
a monopole antenna as the receiver. The PC used for 
the negative refraction experiment has 9 layers in the 
propagation direction (TM) and 33 layers in the lateral 
direction (TM). The interfaces are along the TK direc- 
tion. In all the experimental and theoretical results, the 
electric field is kept parallel to the bars. The operational 
frequency that leads to a structure of approximatively 
n = —1 is 6.5 GHz which correspond to a dimensionless 
frequency of fa/c = 0.345. For the negative refraction 
the horn antenna is oriented such that it makes an angle 
with the normal to the TK interface. 

To examine the negative refraction, we first measure 
the transmission along the first interface in the TK di- 
rection without the PC. We repeat the measurement but 
this time with the PC and calculate the transmission in 
the vicinity of the second interface. The results are plot- 
ted in Fig. 3. It is shown clearly how the center of the 
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Figure 2: ( color online) (a) The photonic band structure 
for the TM polarization, (b) The Equal Frequency Surfaces 
(EFS) with different dimensionless frequencies shown on the 
corresponding curves. 



outgoing Gaussian beam is shifted toward the left-handed 
side of the center of the incident beam. This is a clear 
indication of the negative refraction occurring inside the 
structure. To get a further insight about the negative or 
positive shift, Fig. 4 explains in details in which cases we 
have positive or negative refraction. 

It is important to notice that experimentally we did 
observe the negative refraction for different angles. This 
confirms that the negative refraction relies entirely on the 
LH behavior of the system and not on the convex shape 
of the EFS£i For a particular angle of = 30°, the 
lateral shift was 3.2 cm which correspond to a negative 
angle inside the LH PC of -15°. 

Further insight about the field distribution inside this 
structure can be gained by observing the E field inside 



ffl 

o 



CO 

Eh 



-10 - 



-15 



-20- 



-25- 



-30 





fc 1st interface 




^« "2nd interface 


' / N 




/ / 


\ \ 


/ >/ 


\ X. 


/ f 




/ negative 


positive ^ — 



-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 



Position (x/A) 

Figure 3: The transmission at the first ( solid line) and second 
interface ( dashed line) of the PC versus the lateral position. 
Arrows indicate the refracted beam's direction depending on 
positive or negative shift. 



Schematics of Positive and Negative Refraction 



Negative Refraction 
Second interface 



Positive Refraction 




First interface 



Figure 4: ( color online) Schematics of positive and negative 
refraction. For negative refracted beam, the center is shifted 
to the left hand side of the incident beam and for the positive 
refracted beam, the center is shifted to the right hand side of 
the incident beam. 



a LH PC's slab plotted in Fig. 5 after 77.33T, where 
T = 2tt/uj is the period of the incident waves. The struc- 
ture length along the lateral direction is about 13. 8 A and 
about 2.98A in the propagation direction. A Gaussian 
beam source with A = 46 mm is placed outside the PC 
making an angle 6 = 30° with respect to the normal to 
the interface which is the TM direction. The FDTD sim- 
ulation is performed by discretizing the real space into a 
fine rectangular grid of (a/54 and a/92 for the x and y 
axes respectively). A time step of St = 24.90 ps is used. 
The distribution of the field inside the PC shows how 
the beam is negatively refracted inside the PC and pos- 
itively refracted outside the PC parallel to the incident 
beam. The theoretical result is in a good agreement with 
the experiment since it gives a negative refracted angle 
of about -15°. 

An important advantage of this structure is its prop- 
erty to allow single beam propagation. The structure 



Figure 5: (color online) The E field of a Gaussian beam 
propagating along 30° direction after 77.33T. The structure 
is 13. 8A in the lateral direction and 2.98A in the propagation 
direction. 
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Figure 6: The transmission along TM direction versus the 
dimensionless frequency. The solid arrows indicate the theo- 
retical gaps as calculated from the band structure. 



was designed in such a manner that only one beam prop- 
agation is allowed eliminating the undesired Bragg re- 
flections inside the PC 5 . One of the interesting features 
of this structure compared to the other LHM based on 
proposed PCs is its high transmission. Fig. 6 displays 
the transmission versus the dimensionless frequency. The 
experimental curve shows the high transmission of the 
structure over a wide range of frequencies including the 
operational one. The arrows in this figure show the the- 
oretical indication of the gaps. As predicted by the band 
structure, the experimental curve shows two gaps in the 
TM direction one between 0.25 and 0.276 and the other 
one starting at 0.41. Thus, the agreement between the 
band structure and the experimental curve is very good. 

In order to get an idea about what is happening in- 
side the structure we examine the transmission in both 
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Figure 7: The transmission versus the frequency for both YM 
and YK direction. 



directions. Fig. 7 shows the transmission in the YM and 
YK direction versus the frequency. This structure at the 
dimensionless operational frequency (0.345) has a good 
transmission in the YM direction and a weak transmis- 
sion in the YK direction. A weak transmission in the YK 
direction enhance the transmission in the propagation di- 
rection which will be a critical issue for the superlensing 
phenomenon. It was arguec^^ that the focusing seen 
at the edge of the first band in a square lattice! was not 
due to negative refraction but to anisotropic propagation 
resulting to the funneling effect. This is not the case in 
our structure. As we will show latter, the focusing and 
the superlensing in our structure relies on negative re- 
fraction and on coupling to the surface states. 
In order to investigate the superlensing phenomenon, we 
experimentally measure the intensity of the focus for the 
same structure used in the transmission measurements ( 
33 layers in the lateral direction and 9 layers in the prop- 
agation direction). Two probes are used. The first one 
is placed at a distance of 0.2 A from the first interface of 
the PC and the second one at the same distance from 
the second interface. We first measure the transmission 
with and without the PC . Fig. 8 shows the two resulting 
curves. This figure demonstrates the high intensity of 
the focus and that the transmission is enhanced by more 
than 20 dB for the setup with the PC compared to the 
one without the PC. The solid curve also shows a high 
transmission at the edge of the structure at x/X ~ 1.75. 
This might due to the limited lateral length of the struc- 
ture and to the support station which might give some 
reflection at the edge. 

Theoretically, we investigate the superlensing within a 
slab of PC with a dimension of 20. 7A in the lateral di- 
rection and 1.20 A in the propagation direction. A point 
source with a dimensionless frequency 0.345 placed at a 
distance 0.15 A from the first interface of the PC is used. 
A time step of St = 26.5 ps and a fine rectangular dis- 



Figure 8: The transmission versus the lateral position of a 
probe placed at a distance 0.2 A from the first interface without 
the PC (dashed curve ) and with the PC (solid curve) at 
6.5 GHz or a dimensionless frequency (fa/c) = 0.345. 
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Figure 9: (color online) Snapshots of the E field after 68.49T 
of a point source located at 0.15 A and its image through a 
PC slab of 20. 7A in the lateral direction and 1.20A in the 
propagation direction, (a) shows the complete structure and 
(b) shows the structure with a surface cut of 0.10a. 



cretization mesh of a/40 and a/92 in the x and y direction 
are used. Fig. 9(a) and Fig. 9(b) show the snapshots of 
the E field after 68.49T for the complete structure and 
the one with 0.10a termination. By the termination or 
cut 0.10a we mean that we cut 0.10a from the first and 
last rows of the bars. Notice that in order to excite the 
surface waves at the interface, we studied different ter- 
minations of the surfac e) 16 ! 17 . Among the different cuts, 
one surface termination succeeds in exciting the surface 
waves as shown in Fig. 9b. In this case the cut was about 
0.10a . For both Figures ( Fig. 9a and 9b) a focus is ob- 
tained. However, a better image with a higher intensity 
is obtained when surface waves are involved. The surface 
waves in Fig. 9b are clearly shown as propagating along 
both the YK interfaces. 

One way to verify the good quality of the image com- 
pared to the source is to plot the average field inten- 
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Figure 10: The normalized average field intensity at the 
source (black curve) and image plane (dashed curve) versus 
the lateral direction. 



sity over a period at both the source and image planes. 
Fig. 10 displays the normalized average field intensity 
versus the lateral direction for the cut structure. The 
full width at half maximum (FWHM) of the image beam 
is 0.35A and it is almost the same as the FWHM at the 
source. This result shows that the cut structure, which 
supports surface waves focuses in perfect way the object. 
Thus with a source that is not too close to the structure 
(0.1 5 A) we succeed in exciting the evanescent waves and 
achieving a better focus. 

Apparently our structure with 0.10a cut surface seems 
to be a good candidate for superlensing because we 
achieved a quite high transmission over a wide range 
of angles. Fig. 11 presents the transmission versus the 
angle of incidence for different values of surface termina- 
tions. By trying different terminations we managed to 
optimize the transmission. For the complete structure, 
the transmission start at 80% for normal incidence, it 
gets high for a short range of angles then it drops to less 
than 20% for angles larger than 70°. However the struc- 
ture with a surface termination of 0.10a (see Fig. 11) 
shows the maximum transmission for all incident angles. 
Thus, the transmission is about 95% for angles up to 
40°, and more than 80% for angles up to 65° and more 
than 50% for angles up to 78°. This result shows first 
how the transmission gets enhanced by the means of the 
excitation of the surface waves for a non homogeneous 



structure. Second, high transmission over a wide range 
of incident angles excludes the possibility of the funnel- 
ing effect 14,1 5 and accentuate the fact that the negative 
refraction in our structure is purely a result of negative 
behavior 5 and not a prefered propagation direction. 

In conclusion we have systematically studied a LH 
structure based on a PC with a negative refractive index. 
Our experimental and theoretical results show the nega- 
tive refraction as well as the superlensing phenomenon in 
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Figure 11: Transmission versus the angle of incidence for dif- 
ferent surface terminations. 



this structure. Furthermore, surface waves at the inter- 
face between air and the PC are excited within a specific 
termination of the surface allowing the reconstruction of 
the evanescent waves for better focus. The calculated 
average field intensity in the source plane as well as in 
the image plane shows almost the same full width at half 
maximum demonstrating the perfect image reproduced 
by this new structure. Its high transmission along a wide 
range of angles makes it a good candidate for observing 
LH behavior in PCs. 
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